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We investigate here the Central Limit Theorem of the Increment Ratio Statistic of a 
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The aim of this paper is a simple proof of Central Limit Theorem (CLT in all the sequel) for 
the convergence of Increment Ratio Statistic method (IRS in all the sequel) to a time varying 
Hurst index. 

Hurst index is the main parameter of fractional Brownian motion (fBm in all the sequel), 
Q\ | it belongs to the interval (0, 1) and it will be denote by H in all the following. For fBm, 
the Hurst index drives both path roughness, self-similarity and long memory properties of 
the process. FBm was introduced by Kolmogorov |20] as Gaussian "spirals" in Hilbert space 
and then popularized by Mandelbrot & Van Ness |22| for its relevance in many applications. 
However, during the two last decades, new devices have allowed access to large then huge 
datatsets. This put in light that fBm itself is a theoretical model and that in real life situation 
the Hurst index is, at least, time varying. This model, called multifractional Brownian motion 
(mBm) has been introduced, independently by Levy-Vehel & Peltier [21] and Benassi et al [9]. 
^ ■ Other generalizations of fBm remain possible, for e.g. Gaussian processes with a Hurst index 
depending of the scale, so-called multiscale fBm [5], when H is piecewise constant as in the Step 
Fractional Brownian Motion see [3], or a wide range of Gaussian or non-Gaussian processes 
fitted to applications (see for example |14|. H]). However, for simplicity of the presentation, in 
this work we restrict ourselves to mBm. 

In statistical applications, we estimate the time varying Hurst index through a CLT. Actually, 
CLT provides us confidence intervals. Different statistics can be used to estimate the Hurst 
index. Among the popular methods, let us mention quadratic variations, generalized quadratic 
variations, see [HI [151 H2L and wavelet analysis, see e.g [Tj or [6]. Above methods can be 
expansive in term of time complexity. For this reason, Surgailis et al [27] and Bardet & Surgailis 
[7] have proposed a new statistic named increment ratio which can be used for estimating the 
Hurst index H and is faster than the wavelet or the quadratic variations methods, at the price 
of a slightly larger variance. 

CLT for the different estimators of Hurst index are presently standard in the case of fBm, but 
became very technical in the case of mBm. The main novelty of our work is the simplicity of 
the proofs. In our point of view, mBm is a fBm where the constant Hurst index H has been 
replaced by time varying Hurst index. It is well known that the random field (H, t) i— >■ B{H, t) 
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is irregular with respect to time t, actually with regularity H which belongs to (0,1). It is 
less kown that this field is infinitely differentiable with respect to H, see Meyer-Sellan-Taqqu 
(1999) and Ayache and Taqqu (2005). Thus, for all time to, we can freeze the time varying 
Hurst index, and the mBm behaves approximatively like a ffim. Eventually, CLT for mBm 
follows from CLT for ffim combined with a control of "freezing error". This new and natu- 
ral technology allows us to go further and obtain for example a CLT for the Hurst function 
evaluated at a finite collection of times and also quantitative convergence speed in the CLT. 
Note that, up to our knowledge, the "freezing Hust index" strategy for estimation in mBm 
was introduced, without further proof, in Bertrand et al [lQj. 

The remainder of this paper is organized as follows. In Section [TJ we recall a definition of ffim 
and the definition of the Increment Ration Statistic. Next, in section[2l we review definitions of 
ffim and mBm and precise the localization procedure (or freezing). The main result is stated 
in Section [H] and some numerical simulations are presented in Section 2J All technical proofs 
are postponed in Section 5. 



1 Recall on fBm and Increment Ratio Statistic 

In this section, we present the Increment Ratio Statistic (IRS) method obtained by filtering 
centered Gaussian processes with stationary increments. Before, we recall definition of the 
processes under consideration. 

1.1 Definition of fBm and Gaussian processes with stationary increments 

We describe ffim through its harmonizable representation. However, it is simpler to adopt a 
more general framework and then specify ffim as a particular case. Let X = (X(t),t G [0, 1]) be 
a zero mean Gaussian process with stationary increments, the spectral representation theorem 
(see Cramer and Leadbetter |18| or Yaglom [28]), asserts that the following representation is 
in force 

X(t) = f(l- e**) ■ f/\0 dW(0, for all t G [0, 1], 

JR 

where W(dx) is a Wiener measure with adapted real and imaginary part such that X(t) is real 
valued for all t. The function / is a Borelian even, positive and is called spectral density of X. 
To insure convergence of the stochastic integral, / should satisfies the condition given by 

/ (iAiei 2 )-/(o^<oo. (i) 

JR 

Example: Fractional Brownian motion with Hurst parameter H G (0, 1) and scale parameter 
a > corresponds to a spectral density given by 

/(£) = a 2 |er (2H+1) for all £ G R. (2) 

In this paper, we denote ffim by B(H,t) when a = 1. Stress that this choice is not the 
conventional one. But, IRS is homogeneous and does not depends on a multiplicative factor. 
Thus, in sake of simplicity, we can impose the extra condition a = 1. 

1.2 Definition of the a-Generalized increments 

In all the sequel, we consider the observation of the process X at discrete regularly spaced 
times, that is the observation of (X(to), . . . ,X(t n )) at times = k/n. Secondly, we consider 
a filter denoted by a of length L + 1 and of order p > 1, where p < L are two integers. 
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It corresponds to an arbitrary finite fixed real sequence o := (oo, . . . , ol) £ having p 

vanishing moments, i.e., 



f if t € {0,...,p- 1} 

J^oj^^O if * = p. 



Consequently, it is easy to prove, for any integer m, that 



L L 
l 1= l 2 =0 



( 



if m € {0, . . . , 2p - 1} 



5^ajZ p ) t^O ifm = 2p. 



^ \z=o 



(3) 



(4) 



The family of such filters will be denotes A{p,L). Then, the a-Generalized increments of the 
discrete process (X(tk)) 0<k<n are defined, for all < k < n — L — 1, as follows 



A a X(t k ) = J2 a i x (tk+l) 

k=0 

and their harmonizable representations are given by 

A a X(t k ) = [ e ^a(-£A0/ 1/2 (£) aw(0 

JR 

where g a (-) is specified as follows 

L 

ila 



(5) 



9a{u) := 



(6) 



Examples: In the simple case where a := {clq = l,Oi = — 1), the operator A a corresponds 
to a discrete increment of order 1, and when a := (clq = l,a\ = —2,a>2 = 1), the operator A a 
represents the second order differences. 

1.3 Definition of the Increment Ratio Statistic 

Let (A a X(tfc)) 0<fc<n _ L _ 1 be the a-Generalized increments sequence defined by ([5]) from the 
discrete observation (X(tk)) 0< f :<n - Then, the IRS introduced by Bardet and Surgailis [7] is 
given by 



n-L-l 



IRS 



a<n (X) = —!— V ^(A a X(t k ),A a X(t k+l )) 
n — L ^-^ 



(T) 



fc=0 



where ?/>(•,•) is described as follows 



\x + y\ 



if (x,y) £l 2 \{(0,0)} 



ip(?,y) := ^ \x\ + \y\ 

1 if (x,y) = (0,0). 
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IRS of fractional Brownian motion 



In the case of the fBm with Hurst parameter H G (0,1), i.e X(t) = Bn{t), Bardet and 
Surgailis have established in [7J Corollary 4.3, p. 13], under some semi-parametric assumptions, 
the following CLT for the statistics IRS ajn 

v^(IRS a ,„(i^)-A a (iT)) 4 A/-(0, £*(#)) W ith | H H l { ^ ] llZ^Z^!) ( § ) 
v 

where the sign — > means convergence in distribution, 



A a (H) := A (p a (H)) (9) 
Ao(r) := ±arc<x)s(-r) + ~^^log(^) (10) 
f 2 2H " 1 -l if a = (1,-1) 

Pa(H) = _ 3 2H + 2 2H + 2 _ 7 

( 8 _ 2 2g+i if a = (1,-2,1) 

and the asymptotic variance T,^(H) is given by 

E 2 a (H) =J2 co v(i>(KBH(t ),A a BH(t 1 )),^(A a B H (t j ),A a BH(t j +i))). 

The graphs of Aq(p), p a (H) and A a (H), with a = (1,-1) or a = (1,-2,1), are given in Fig- 
ure [H Figure [2] and Figure [3] below. It is easy to prove that the function H i— > A a (H), with 



Figure 1: The graph of Ao(p). 



Figure 2: The graphs of p a {H) with a = (—1, l)(left) and a = (1, —2, 1) (right). 



a = (1, —2, 1), is a monotonic increasing function in the interval (0,1), see Figure[3j Therefore, 
H n = A" 1 (IRS a . n (-Bf/)) provides an estimator of the Hurst parameter H with convergence 
rate O {\fn). Moreover, we refer to Stoncelis and Vaiciulis |26| for a numerical approximation 
of the variance Y!^(H) with a = (1, —1) or a = (1, —2, 1), needed for construction of confidence 
intervals, see [7J Corollary 4.3, p. 13 and Appendix, p. 32]. 
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Figure 3: The graphs of A a (H) with a = (— 1, l)(left) and a = (1, —2, 1) (right). 



2 Going from fBm to mBm and return by freezing 

The main goal of this section is to present different representations for the FBm and the mBm 
enabling us to present the time freezing strategy we will use to prove our main theorems. 

FBm and its different representations 

Fractional Brownian motion was introduced by Kolmogorov |20| and then made popular by 
Mandelbrot & Van Ness [22j. This process has been widely used in applications to model data 
that exhibit self-similarity, stationarity of increments, and long range dependence. FBm with 
Hurst parameter H £ (0, 1), denoted by (Bi{(t),t £ [0, 1]), is a centered Gaussian process with 
covariance function defined for s, t £ [0, 1] by 

E [B H (t)B H {s)] = \ {t 2H + s 2H + \t- s\ 2H ) . (12) 

This process is characterized by its Hurst index which drives both pathwise regularity, self- 
similarity and long memory, see e.g. the overview in Bertrand et al (2010). Before going 
further, let us precise notations: in all the sequel we will denote by Bh the fBm and B(H,t) 
the random field depending on both Hurst index and time. Up to a multiplicative constant the 
two notions coincide, more precisely we have Bn(t) = C(H) x B(H,t) for C(H) a non-negative 
constant depending on H. 

Fractional Brownian motion, (Bjj(t),t £ [0,1]), can be represented through its harmonizable 
representation (JTJ [2j) , or its moving-average representation (see Samorodnitsky & Taqqu |25l 
Chapter 14]). A third representation is the wavelet series expansion introduced by Meyer et al 
[23], and then nicely used by Ayache and Taqqu (2003 and 2005). In this case, it is judicious 
to shift to the random field representation defined as follows 

B(H,t) = J2J2 a ^ H ^> e ^ for all i € [0, 1] (13) 
jez fcez 

where (£jk)(j,k)ez 2 is a sequence of standard Gaussian random variables A/"(0, 1), the non- 
random coefficients ajk(t,H) are given by djk(t,H) = I (1 — e'*^) • |CI~^ +1 ^ 2 ^V'i,fc(0 d£, and 

t/jj^k is the Fourier transform of the Lemarie-Meyer wavelet basis ipj k- Let us refer to Ayache 
and Taqqu (2003) for all the technical details. To put it into a nutshell, by using the Meyer 
et al 's Lemma (|23|). we can prove the existence of an almost sure event Q* , that is such 
that Pr(0*) = 1, such that for all u £ fi* the series (fl"3l ) converges uniformly for (H,t) £ K 
where K is any compact subset of (0, 1) x R. Moreover, the field is infinitively differentiable 
with respect to H with derivatives bounded uniformly on every compact subset of (0, 1) x R 
by a constant C*(oj) > where C* is a positive random variable with finite moments of every 
order. 
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MBm and its different representations 

To be short, mBm is obtained by plugging a time varying Hurst index t t— > H(t) into one of 
the three representations of the fBm given above, that is the moving average representation, 
the harmonizable one ((UJ [2])) or the wavelet series expansion (|13l ), The function t i— > H(t) 
should be at least continuous, and if the Holder regularity of the function t i— > H(t) is greater 
than max(H(t),t 6 [0, 1]) (the so-called condition (C) in Ayache and Taqqu (2003)), then for 
every time t G [0, 1] the roughness of mBm is given by H(t). Les us also refer to Cohen |17| 
where he proves that the moving average representation and the harmonizable representation 
of mBm are equivalent up to a multiplicative deterministic function, and to Meyer et al. to 
the almost sure equality of harmonizable representation and wavelet series expansion. 

With this tools, we are now in order to precise our "freezing" technology: 



MBm behaves locally as a fBm 

By applying Taylor expansion of order 1 around any fixed time t* S [0,1], we obtain the 
following formula 

B(H(t),t)I n * = B(H*,t)l u *+R(t)l u * (14) 

where R(t) refers to the Taylor rest which satisfies 

sup |i2(s)I Q *| < C*(ui)\H(t) - H*\ (15) 
se[o,i] 

with C* > a positive random variable with finite moments of every order. Noting that H* 
corresponds to the value of the Hurst function H(-) at t* , i.e., H* = H(t*) and Iq* represents 

the indicator function of a subset $7* defined by: Iq* (uj) = < ^ ^ w ^ ^ _ Next, if we know 



else 

that the Hurst function H (•) has a Holder regularity of order i] > 0, so we obtain immediately 
that 

sup \R(s)I n *\ <M*I n *|i-i*|" (16) 
se[o,i] 

with M* = c x C* > a positive random variable with finite moments of every order. 



3 Main results 

This section is dedicated to the CLT of the IRS localized version for the mBm. Let us how- 
ever first give a simple result on the CLT for the IRS of Gaussian processes with stationary 
increments, which, applied to the fractional Brownian motion, gives with a simple proof the 
result of Bardet-Surgailis [7]- 

We thus consider a process X observed through the knowledge of (X (to) , . . . , X (t n )) with 
tk = k/n for k = 1, . . . ,n. The corresponding increment ratio statistic IRS ajn (X) is defined 
by ©, with a filter a G A(p, L), that is satisfying ([3]). 

Theorem 1 (Fractional Brownian motion) 

i) Let X be a zero mean Gaussian process with stationary increments. We assume that 

EMi)l 2 <+°° ( 17 ) 

where r a (j) := cov(A a X(to), A a X(tj)) for j G Z is supposed independent of n. Then 
^(IRS a , n (X)-Ao( Pa )) 4 N(0,Z 2 a ) (18) 
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where Ao(-) is defined by ()10p . p a represents the correlation between two successive a- 
Generalized increments, and the asymptotic variance is given by 

E 2 a := £ cov(V (A a X(t ), A a X(ti)) , $ (A a X(t,), A a X(i, +1 ))) 

and is weJI defined and belongs to [0, +oo). 

ii) In particular, let X be a fBm, that is X = Bu with Hurst parameter H G (0,1). 
Moreover, in the case p = 1 assume the extra assumption H G (0, 3/4). Then CLT U8\) 
is in force where A a (-) = Ao(/9 a ( - )) JS a monotonic increasing function of H with Ao(-), 
resp. p a { ) described by (|10|) . resp. (|25p . and the asymptotic variance T, 2 is given by 

Y 2 := cov(^ (A a B H (t ), A^fa)) , ^ (A a fl H fe), A a B# (i i+ i))) 
jez 

which is well defined and belongs to [0, +oo). □ 
Remark: 

1. In the sequel, in order to could inverse function A a (-), we will assume that the filter 
a G A(p, L) satisfy L = p > 1 and a; = (— l) p ~ l (^) for all I e {0, . . . ,p}. The class of such 
filters will be denoted B{p) and named binomial filters. This restriction is motivated by 
the fact that in the particular case where a G B(p), the correlation function p a (-) defined 
by (|25p . is a monotonic increasing function of H, instead of in the general case where 
a G L) it is not always true. 

2. The regularity of A a (H) enables us then to get via the well known Delta- method the 
CLT for the Hurst parameter. However no closed formulae for A a (if)~ 1 is available so 
that the limiting covariance will be no further explicit. 

3. We stress once again that the proof of the theorem is quite simple. Note also that using 
recent results of Nourdin et al |244 Th. 2.2], we even have that there exists a sequence 
7(n) decaying to zero such that for all h G C 2 and ~ A/"(0, 

\E[h(y/^(lRS a , n (X)-A (p a ))) - MAO] | < Halloo 7(")- 

The precise estimation of j(n) is however out of the scope of the present paper and will 
be found in |19| . Using |24[ Cor. 2.4], we also have that the previous CLT may be 
reinforced to a convergence in 1-Wasserstein distance or in Kolmogorov distance. 

4. The reader will have noticed that the assumption r a (j) := cov (A a X(to), A a X(tj)) in- 
dependent of n for all j G Z is a quite strong one. Indeed, for the multiscale Brownian 
motion, this not true. However, in a sense, it is asymptotically true and it may then 
be applied to prove the convergence of the IRS to the Hurst parameter related to the 
highest frequency. See [19] for further details. 

To achieve our final goal, we state by presenting a Lemma where we prove that the localized 
version of IRS for mBm converges in L 2 (0) to the IRS of fBm with a certain rate. 

Localized version of the IRS for mult ifr actional Brownian motion 

Let us consider a multifractional Brownian motion with Hurst function H(-) denoted by X = 
(B(H(t),t),t G [0, 1]). Secondly, let t* G (0, 1) be an arbitrary fixed point, then we denote by 
v n (lit*) the set of indices around t* , given by 



i/ n (7,0 = {ke{0,...,n-L-l}:\t k -t*\<n-i} 
= {[nt* -n 1 ~"<\,...,\nt* + n 1_7 J} 



(19) 
(20) 
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where \_x\ is the integer part of x and 7 £ (0, 1) is a fixed parameter which allows to control 
the size of v n (7, t*) which cardinal is equal to v n {^) := 2n 1_7 + 1. Finally, for any n large 
enough, we denote by IRS^ (B H ^ the localized version of IRS defined as follows 

IHS# (**(•)) = 2w i- 7 + 1 E ^(A^^ftO.AaB^,^!)). (21) 

fc=[ni*-n 1 -fj 

With these notations, we are in order to state our main result: 

Theorem 2 ( Mult ifr actional Brownian motion) 

i) Let X be a mBm, its localized IRS be defined by (121)) and assume that 7(1 + 77) > 1. 
Then 

n d-7)/a ( ffiS 7j (BffH ) _ Aa(ir) ) 4 aa( ,^) with { f;^^( 4) ;^>2< 22 ) 

where A a (-) = Ao(/3 a ( - )) JS a monotonia increasing function of H* with Ao(-) & Pa( - ) 
described by (|10p & (|25p , and the asymptotic variance is given by 

:= £ cov (A a Bn* (to), A n B ff . (t x )) , ^ (A a B H * fa), A a B H * (t j+l ))) 

jez 

and is well defined and belongs to [0, +00). 

ii) Let now consider < to < ti < ... < t m for a finite m then under the same assumption 
we can enhance the previous CUT to the vector 

n^' 2 (lRSy$ (-Bff(.)) - A a (H(ti)), IRSlf™ (Bhq) - A a (H(t m ))) 

with a well defined limiting covariance (3. □ 
Remark: 

1. Here again, one can use results of |24j to get explicit estimates on the speed of convergence 
for this CLT. 

2. It is highly interesting to upgrade the previous CLT to the trajectory level, needing then 
a tightness result, for example to test if the Hurst coefficient is always greater than 1/2, 
or to perform other test. Such kind of result will be developed in |19] , 

4 Numerical results 

In this section, for numerical estimation of the Hurst index by IRS, one has chosen a binomial 
filter of order 2, i.e. o = (1,-2,1), insuring the convergence of the estimator H n for any 
H € (0, 1). At first, we analyze through Monte-Carlo simulations the efficiency of the Hurst 
parameter of fBm estimator given by IRS. Then, we study the estimators of some Hurst 
functions of mBm obtained by localized version of IRS, and we compare it with the estimators 
given by Generalized Quadratic Variations (GQV) method, see e.g Coeurjolly |16| . 

Estimation of the Hurst index of fBm 

At first, by using Wood and Chan |13] algorithm, for n = 10000 we have simulated three repli- 
cations of the fBm sequences Bh = (Bfj(to), . . . , i?//(t n )), at regularly spaced times such that 
tfc = k/n with k = 0, . . . , n , for three values of the Hurst parameter H, denoted {Hi, H2, H3}, 
and given by 
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(Ci_) Hi = 0.3 < 1/2 for short range dependent case, 
(C2) H2 = 1/2 for standard Brownian motion, 
(C3) H3 = 0.7 > 1/2 for long range dependent case, 
see Figure H] below. 

Then, for each sample (C;) with i € [1,3], we have computed the increment ratio statistic 




Figure 4: Simulated fBm process with H = 0.3 (left), H = 0.5 (middle) and H = 0.7 (right). 

IRS2. n (BHi) an d estimated the Hurst index given by H n ^ = A^" 1 (IRS2,n(-Bf/i))- We remark 
that the IRS methods provide good results given in Table [1] below. 
These examples are plainly confirmed by Monte Carlo simulations. Indeed, for each case [C\) 



Exact values of H 


0.3 


0.5 


0.7 


Estimated values of H 


0.3009 


0.4993 


0.7000 



Table 1: Estimated values of H . 



with i G [1,3], we have made M = 1000 simulations of independent copies of fBm sequences 
= (B^(to), B$(tnj) , for k = 1, . . . , M. We find also good results illustrated by 

the following histograms, see Figure [5j which represent the distribution of the estimator H n ^, 
for % € [1, 3]. Thus, we have computed the standard deviation E|iT— H\ 2 given in Table [2] below. 




Figure 5: Distribution of the estimated values of H in the case H = 0.3 (left), H = 0.5 (middle) 
and H = 0.7 (right). 
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Different values of H 


0.3 


0.5 


0.7 


Standard deviation 


8.1865 x 10~ 5 


7.6914 x 10~ 5 


6.9837 x 10" 5 



Table 2: Standard deviation. 



Local estimation of the Hurst function of mBm 

To synthesis a sample path of a mBm, one has used the Wood and Chan circulant matrix 
improved with kriging interpolation method, which is faster than Cholesky-Levinson factoriza- 
tion algorithm. In fact, both methods are not exact but provide good results. For, n = 10000, 
we have simulated three samples of the mBm sequences Bh{-) = (-^i/(t )(*o)) • • • j Bjj(t n ) (tn)) j 
at regularly spaced times such that i& = k/n with k = 0, . . . , n , for three types of the Hurst 
function H(-), namely 

(C4) Linear function: H±(t) = 0.1 + 0.8i, 

(C5) Periodic function H^(t) = 0.5 + 0.3sin(7ri), 

0.3 

(Cq) Logistic function: Hq{€) = 0.3 + 



(l + exp(-100(t-0.7)))' 
see Figure [6] below. 

Then, for each sample (Ci) with i £ [4, 6], we have estimated the Hurst function H n { by using 



Figure 6: Simulated mBm process with H(-) linear function (left), H(-) periodic function (mid- 
dle) and H(-) logistic function (right). 



the localized version of IRS with 7 = 0.3 and the GQV method. We note that both methods 
provide correct results represented by Figure [7] below. 

These results are plainly confirmed by Monte Carlo simulations. Actually, for each case (Ci) 
with i G [4,6], we have made M = 1000 simulations of independent copies of mBm sequences 

Bff\\ = (B^^(to), . . . , Bj£}, t ^{t n )J , for k = 1, . . . , M. Then we have computed the Mean 



Integrate Square Error (MISE) defined as MISE = E \H(tj) - H(t 



\ 3=0 



E\\H- 



-^1Il 2 (o 1) wn i c h i s a criterion widely used in functional estimation, see Table [3] below. 
We observe through Table [3] that both methods provide globally the same results when the 
function H(-) varies slowly (see linear and periodic cases), whereas in the case where H(-) 
presents the abrupt variation it appears that the GQV is a bit more precise compared to the 
IRS method. 
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Figure 7: Estimation of the Hurst function H(-) with H(-) linear function (left), H(-) periodic 
function (middle) and H(-) logistic function (right). The graphs of function H(-), Hn (•) 
its estimation by IRS, and H^® V \-) its estimation by GQV, are represented in green, red and 
blue respectively. 





iJ(-) Linear 


H(-) Periodic 


H(-) Logistic 


MISE by IRS 


2.6743 x 10" 4 


1.4743 x 10" 4 


5.3546 x 10~ 3 


MISE by GQV 


8.9547 x 10~ 4 


5.4743 x 10~ 4 


8.9743 x 1(T 4 



Table 3: MISE given by IRS method and GQV method 



5 Proofs of the main results 

This section contains the proof of the results of Section 3. Note that we have divided the 
proof of Theorem 1 in two parts: first we consider the general case of Gaussian processes with 
stationary increments and then in a second part we investigate the application to fractional 
Bronwnian motion. 

5.1 Proof of localization 

First, we can deduce as a corollary that the a-Generalized increments sequence (A a X(tfc)) 0<fc<n _ i _ 1 
form a family of stationary identically distributed centered Gaussian r.v. with variance 

al n = cov(A a X{t k ),A a X(t k )) 
5a(e/n)| 2 -/(0 d£ 

\9a(Un)\ 2 -f(0dC, 

covariance given, for all < k±, ki < n — L — 1, by 

cov a) „(t fel ,t fc2 ) = cov (A a X(t kl ), A a X(t k2 )) 

^ kl - k2 ^ /n \9a^/n)\ 2 ■ f(0 d£ 



cosa^-W^ls^/n)! 2 -/^) dd, 

and correlation between two successive a-Generalized increments defined by 

cov (A a X(t k+1 ), A a X(t k )) cov a , n (*fe+i, t k ) 



Pa,r 



[cov (A a X(t k+1 ), ApXit^))} 1 ' 2 ■ [cov (A a X(t k ), A p X(t k ))] 1/2 <n 



Pierre, R. BERTRAND, Mehdi FHIMA and Arnaud GUILLIN 



12 



where g a (-) is described by ([6]). Therefore, for a fixed < k < n — L — 1, it is easy to remark 
that there exist two independent standard Gaussian r.v. Z k ,Z k+ \ ~ AA(0, 1) such that 

A a X(t k ) = a a>n Z k (23) 



A a X(t k+1 ) = a a , n (p atn Z k + Jl-pl n Z k+1 ), (24) 



where the sign ~ means equal in distribution. 

Remark: In the particular case of the fBm, the correlation between two successive a-Generalized 
increments, p a , n , does not depend on n. Indeed, we know that the spectral density of the fBm 
is given by Q, then we have 

/R + cos(Wk(W| 2 -r (2g+1) dj 

PaA) ; R+ k(e/n)i 2 -e-^ +1 )de ' 

And after, we can change variable ^/ntoti. So this implies that 

f R+ cos(u)\g a (u)\ 2 • n-^+D du _ 

which is independent of n. 

5.2 Proof of CLT for Gaussian Processes with stationary increments 

The proof uses the notion of Hermite rank and Breuer-Major theorem, see for e.g Arcones [21 
Theorem 4, p. 2256] or Nourdin et al [24| Theorem 1, p. 2]. 

Definition 1 (Hermite rank) 

Let G be a R rf Gaussian vector and (ft : R rf — > R be a measurable function such that E |(/>(G)| 2 < 
+oo. Then, the function (f> is said to have Hermite rank equal to the integer q > 1 with 
respect to Gaussian vector G, if (a) E [(<j)(G) — E (<j)(G))) P m (G)] = for every polynomial P m 
(on R d ) of degree m < q — 1; and (b) there exists a polynomial P q (on R rf ) of degree q such 
that E [(0(G) - E (^(G))) P q {G)} / 0. 

We first give the proof of Theorem 1 in the general case of Gaussian processes with stationary 
increments and then in a separate part the application to fractional Brownian motion. 

Proof of Theorem 1 

First, in the sequel we denote by G = (G k , G k +\) = (A a X(t k ), A a X(t k +±)) these two successive 
stationary a-Generalized increments defined by ( ([23p . (|24[) ) . Then, according to Bardet and 
Surgailis Appendix, p. 31], we know that 

E[iP(G k ,G k+1 )] = A (Pa), and E |^ (G k , G k+l )\ 2 < +oo, 

where Ao(.) is defined by (fTU|) and p a is the correlation between G k and Gfc+i. To achieve our 
goal, we start by defining a new function cp : R 2 — > R such that 

( p(X,Y)=i>(X,Y)-A (p a ). 

Then, (J) is in fact a Hermite function with respect to Gaussian vector G = (G k ,G k +\) with 
rank equal to 2. Therefore, by applying Breuer-Major theorem, see e.g Arcones [21 Theorem 
4, p. 2256] or Nourdin et al [241 Theorem 1, p. 2], we get directly the CLT (|18[) . So, the key 
argument of our proof is to determine the Hermite rank of (p. We include here the proof of 
the fact that the Hermite rank is 2 as the proof does not seem to appear elsewhere. Let 
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Pq(X,Y) = c ,P 1 (X,Y) = c u X + ci 2 Y + cio and P 2 (X,Y) = X 2 be three polynomials (on 
M?) with degree respectively 0, 1 and 2. First, it is easy to see that E [4>(G)Pq(G)] = 0. Now, 
we must to show that E y>(G)P 1 (G)] = 0. We have 

E [0(G)Pi(G)] = c n E [4>(G)G k ] + c 12 E [0(G)G? fc+ i] + c 10 E [0(G)] . 

=0 

Then, 

E [0(G)P a (G)] = c n E [^(G fc , G fc +i)G fc ] - en A (p a ) E[G fc ] +c 12 E [0(G fc , G fe +i)G fe +i] - c 12 A (p a ) E[G fc _ 



because G^ and G k+ \ are zero-mean r.v. and due to the fact that the r.v G k and G k+ \ have 
a symmetric function, we can write without any restrictions that 

E[0(G)P(G)] = (cn + c 12 )E [0(G fc ,G fc+1 )G fc ]. 

By using definition of (G k ,G k+1 ) = (A a X(t k ), A a X(t k+1 )) given by (([23]), ®), we get 



(en +ci 2 )- 1 E[0(G)Pi(G)] = a a E a a Z fc ,a a (p a Z k + ^1 - p 2 Z k+1 )) Z k 



where Z k and Z k+ \ are two independent standard Gaussian r.v. Z k ,Z k+ i ~ M (0,1). Thus, 
by using homogeneity property of •) specified by: ijj(aX, aY) = ip(X,Y), we obtain 



(en + c^)- 1 ! [0(G)Pi(G)] = a a E ^ (z fc , (p a Z fc + Vl-pS^fc+ij) Z fc 
Next, we have 



(c 11 + c 12 )^ 1 E[</)(G)P 1 (G)] 



2tt 



-0 [Zl,p a Zi + ~ Pi z 2 ) 2l 
\ / 



exp 



z\ + z\ 



dz\dz2 



And after, we can change variables (zi,z 2 ) to polar coordinates (r cos((9), r sin(0)). So, this 
implies that 



2vra- 1 ( Cl i + c 12 )- 1 E [0(G)Pi(G)] 



(cos(0),p o cos(0) + V 1 - Pi sin(0)j cos(9)d9 



r 2 exp 



Then, we remark that 



^ (cos(0),p a cos(0) + a/1 - ^ sin(0)) cos(0)d0 

cos{9)d9 



cos(9) + p a cos(0) + \/l - 


- Pa sm(6») 


|cos(0)| + 


p a cos(0) + a/1 


- Pi sin(0) 



COS (li — 7r) + ,0 a COS (it — 


tt) + a/1 - 


- Pi sin(n - tt) 


COs(li — 7r) 1 + 


p a cos(u - 


-vr) + A/T 


- Pi sm(-u - tt) 



cos(u — ir)du 



cos(it) + /> a cos(u) + a/1 - pi sin(u) 



I cos (it) I + p a cos(u) + a/1 - pi sin(it) 



cos(u)du 



i/> (cos(0),p a cos(0) + ^/l - p2 „ S in(0)) cos(0)d0. 
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Thus, we deduce directly that E [(j){G)P\{G)] = 0. In the similar way, it is easy to prove that 
E [(f>(G)P2(G)] ^ 0. So, by using Definitional we can sa Y that is a Hermite function with 
rank equal to 2. Therefore, Theorem [T] becomes an application of Breuer-Major theorem which 
gives directly the proof of CLT (|18p . ■ 

Proof of CLT for FBm 

Next, we present the correlation function properties of the a-Generalized increments sequence 
of a fBm. 

Property 1 (Correlation function of the a-Generalized increments) 

Let (Bf{(t),t € [0, 1]) be a fBm with Hurst parameter H £ (0, 1) and let (A a Bjj (ife))o<fc<n-£-l 
its a-Generalized increments sequence defined by ©, with a 6 A(p,L) a filter given by ([3]). 
Then, for all j £ 7L, we have 

r a ,n(j) = - 1 2H x C a (j) 
where r a ^ n (j) := cov (A a B H (t ), A a B H (tj)) and C a (j) is given by 



Ca(j) 



^ a h a h \j + l 2 -h 
h,h=o 



\2H 



2H 

2p 



x J 



2H - 2 ? with 



J=0 



2H\ Ufo^H-k) 



2p 



{2p)\ 



And the correlation between two successive a-Generalized increments, is specified by 

Pa(H) 



Ei 1 ,z 2 =o a h a h \l 2 - h\ 2H 



(25) 
□ 



Proof of Property 1 

To compute the covariance function of the a-Generalized increments sequence, we start by 
using the initial formula of the covariance function of a fBm defined by (fT2"|) . Then, we obtain 

r a ,n(j) ■= cov (A a B H (to),A a B H (tj)) 

L 

= ^2 a h ai 2 E[B H (t h )B H {tj+ h )] 
h,h=o 

L 



1 ^ > 2H \ ^ 2H \ ^ i 



211 



h,h=0 



1 



h,h=0 



h,h=0 



=0 



X C a (j) 



2-n 2H 

where C a (j) = l2=0 a^a^j + k — h\ 2H ■ Now, we give an equivalent of C a (j) when j — > +oo. 
To do this, we use the Taylor expansion as follows 

L 



C a {j) = j 2H x a hah 

h,h=0 



1 + 



211 



J 2H x 



L 



Yl a h a h + £ ( 2 f )j k a hai 2 \h-h\ k 

h,l 2 =0 k=l ^ ' iiJ2=0 
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Next, by using (JU), we know that when we sum over k, every term in the expansion gives a 
zero contribution for any integer k < 2p. So this implies that 

This finishes the proof of Property [TJ ■ 

And after, we note that the function satisfies the homogeneity property specified by: 

ifj(aX, aY) = i(j(X,Y). So, this allows us to rewrite IRS ajn (-Bf/) as follows 

1 n— L— 1 

IRS a , n (B H ) = V */j(A a B§ d (t k ),A a B% d (t k+1 )) 

n — L L — ' 

k=0 

where A a Bfj- d represents the standardized version of A a Bu described, for all < /c < n — L— 1, 



as 



a Rstd/, ^ A a B H (t k ) 



y/Var [A a B H (t k )} 
and its covariance function is given by 

r a (j) := cov(A aJ Bg d (t ),A aJ B|* d (^)) 
r a , n (j) 
r a ,nU) 

Ca(j) 



C a (0) 



2 

i /2ir 



i^+oo C a (0) V2p, 

which is independent of n. So, according to Theorem [TJ the key argument is to prove that 

^2K(j)\ 2 < +oo. 

Thus, by using a Riemman sum argument, we can deduce immediately that this condition is 
verified if and only if AH — 4p < —1, i.e H < p — 1/4, and this implies that H G (0,3/4) if 
p = 1 and that H G (0,1) if p = 2. Therefore, the assumption (|17|) of Theorem Q] is satisfied 
and so we obtain a simple intuitive proof of the CLT (|18p applied to the IRS of fBm. This 
finishes the proof of Theorem [TJ ■ 
End of The proof of Theorem 1 



5.3 Proof of CLT for mBm 

The proof of Theorem 2 relies on a localization argument given in the following Lemma 
Lemma 1 

First, we consider t* G (0, 1) be an arbitrary fixed point, 7 (0, 1) a fixed parameter which 

allows to control the size of the indices set around t* , and a G B{p) a binomial filter. 

Let B h q = {B H(t) (t),t G [0, 1]) be a mBm with Hurst function H(-) G C v ([0, 1], [H , H°]) and 
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B H * = (B H *(t),t £ [0, 1]) a fBm with Hurst index H* = H(t*). 

Moreover, we consider IRS^ n the localized version of IRS for mBm defined by ([21 

and IRS atn (Bh* ) a modified version of the IRS for fBm described as follows 

lnt*+n 1 —'\ 

IRS a>n {B H * ) = 2ral „ 7 - 1 * ^B H * (**)• A " Bh * • 

k=[nt* -ni-T J 



Then 

E 



ms£(B H n)-mSa, n (B H *) 2 = o (rr™). (26) 



Proof of Lemma 1 

For n large enough, we have 

2 



I airs:™ i a 



IRS2;* ( J B H (.))-iRS a ,n( J BH*: 



2n 1 -T + 1 



[nt*+n 1 -Tj 

^(A a 5 H(tfc) (i fe ),A a S H(tfc+l) (t fc+ i)) -4)(A a B H .(t k ),A a B H *(t k+1 )) 

fc=Lnt*-n 1 -Tj 



Then, by using Cauchy-Schwarz inequality, we get 

+ ' 



|AIRS™| 2 < - i £ E|^(A o J5 H(tfc) (t0 ) A o J3 H (t w )(*fc+i)) -^(AaB H .(t fc ) J A a J3 H .(t fc+1 ))^ 



AIRSrf < 2nl „ 1 7 + 1 £ |^(A aJ B H(tfc) (t fc ),A Q ^ (ifc+l) (t fc+1 )) -^(A a B H ,{t k ),A a B H *{t k+l ))\ 2 . 

k=[nt*-n 1 -~t] 

This implies that, 

Lnt*+n 1-7 J 

E| 

fc=[nt*-n 1 -Tj 

Now, we recall that f2* represents the event with probability 1 introduced in Subsection 2.3. 
Then, according to Bruzaite & Vaiciulis |12|. Lemma 1, formula 3.3, p. 262] and our Taylor 
expansion (|14|) . we deduce that 

E{|V(A a B H(tjb) (t Jfc ) > A a S ff(tfc+l) (t fc+1 )) -il>(A a B H .(t k ),A a B H .(t k+1 ))\ 2 I a *} 

< C ( Pa (H*)) -E{(\A a R(t k )\ + \A a R(t k+1 )\)I Q *} • 

where the constant C (p a (H*)) > depend only on p a (H*) and (A a R(t k )) corresponds to 
a-Generalized increments at t k of the rest (R(t),t E [0,1]) defined by (|15|) . Next, by using 
([161) . we deduce that there exist a constant k > such as 

E|V(A a fl H(tfc) (t Jfc ),A a B ff(tjk+l) (^ +1 )) -V(A a S ff .(t fc ),A a B ff .(^ + i))| 2 ]I . < ktT™ 
Therefore, we obtain 



Moreover, we know that 
E 



E{|AIRS™| 2 I n *} < wr™. 
AIRS™ | 2 = E { | AIRS™ | 2 1^ } + E { | AIRS™ | 2 In\ n . } , 



=o 

by applying Cauchy-Schwarz inequality. Then, we get directly 



EIAIRS™! 2 < Kn~^. 



This finishes the proof of Lemma [TJ 
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Proof of Theorem 2 

First, according to our Lemma [1] we have 



IRS^(^ ( .))-IRS a , re (^) = O (n"™) 

n— >+oo 



Next, it is easy to see that 



E 



n^' 2 (iRSlf; {B^) - A a (H* 



n 



(l-7)/2 



(IRS , n (B H *) - K(H*)) 



O (V _t(1+ ^ 



After, by applying Theorem [TJ we know that 



n 



(i-7)/2 



(IRS a , n (B H .) -A a (H*)) 4 AA(0,S^) with 



F* G (0, 3/4) if p = 1 
H* e (0, 1) if p > 2 



Therefore, CLT (|22p is satisfied if and only if 7(1 + rf) > 1. 

Let us now sketch how to extend it to the multidimensional case: first we may operate a 
multidimensional freezing of time in the sense that there exists an almost sure event Q* such 
that 

Vt, B(H(t),t)l n * = B\H{t u t)l n * + iT(t)l n *, sup |if(s)I n *| < C*(ui)\H(t) - H*\ 

s6[0,l] 

and the process -B* are defined using wavelet expansion so that the correlations between them 
are well described. We may then consider fractional Brownian motions rather than mBm. 
Secondly we use Cramer- Wold device (see e.g. Th. 7.7 in Billingsley |11| ) : it is sufficient to 
get the CLT for every real numbers b\ , . . . , b m for 



(l- 7 )/2 



i=l 



J> (iRSl^Bij^) - A a (H(U)) 



which is obtained exactly as before. 
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